Extended Hubbard Model with Unconventional Pairing in Two Dimensions 



in 
o 
o 

(N 



Akinori Tanaka^ and Masanori Yamanaka^ 
^Department of Applied Quantum Physics, Kyushu University, Fukuoka 812-8581, Japan 
^'Department of Physics, College of Science and Technology, Nihon University, Tokyo 101-8308, Japan 

(Dated: March 31, 2005) 

We rigorously prove that an extended Hubbard model with attraction in two dimensions has an 
unconventional pairing ground state for any electron filling. The anisotropic spin-0 or anisotropic 
spin-1 pairing symmetry is realized, depending on a phase parameter characterizing the type of 
local attractive interactions. In both cases the ground state is unique. It is also shown that in a 
special case, where there are no electron hopping terms, the ground state has Ising-type Neel order 
at half-filling, when on-site repulsion is furthermore added. Physical applications are mentioned. 

PACS numbers: 71.10.Fd, 74.20.-z 
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Unconventional superconductivity with gap symme- 
tries other than the conventional s-wave has been found 
ubiquitously in correlated electron systems. Examples 
include heavy fermions high Tc cuprates 0, ruthen- 
ate 0, organic conductors etc. The common feature 
of those compounds is proximity of antiferromagnetic or 
ferromagnetic order. A vast number of studies have been 
extensively devoted to revealing nature of these phenom- 
ena. So far, within the mean-field approach, it is rec- 
ognized that effective electron-pair attraction depending 
on electron momentum can cause unconventional super- 
conductivity, but, there still is no convincing evidence 
which model captures the mechanism truly. As a many- 
body problem, it is an extremely hard task to make a 
definite criterion to distinguish the validity of the mod- 
els beyond the mean-field level. Indeed, electron systems 
exhibit various physical phenomena, relying on a subtle 
interplay between kinetic and interaction energies. It is 
thus desirable to rigorously establish the occurrence of 
unconventional pairing in concrete models of correlated 
electrons. The attempts in this direction will shed light 
on the mechanism for unconventional superconductivity 
and, in turn, give us useful information about possible 
sources of effective pair attraction in real materials. 

In this paper, we rigorously construct a series of the 
electronic models having ground states with unconven- 
tional pairing symmetries. We consider a two dimen- 
sional tight-binding model with attractive interactions 
which act on electrons occupying certain localized single- 
electron states (corresponding to and © below) . For 
each even number of electrons, the model is proved to 
have the unique ground state [3 in which all electrons 
form anisotropic pairs with spin or spin 1 , depending on 
a phase parameter (0 in (|3J)) of the localized states. Here, 
we treat the model in two dimensions, since the case is 
most relevant to the experiments mentioned above. Ex- 
tensions of the present method and idea to higher dimen- 
sional systems or other lattice structures are straightfor- 
ward. 

Remarkably, unlike usual mean- field Hamiltonians, our 
Hamiltonian conserves the electron number. The occur- 
rence of the electron-pair condensation is thus non-trivial 
in the present model. To our best knowledge, this is the 



first time that a model Hamiltonian is proved to exhibit 
condensation of unconventional electron pairs including 
spin-1 pairs 0. Further advantage of our model is a re- 
lation to proximity of magnetic orders. In a case, where 
there are no electron hopping terms, the model exhibits 
antiferromagnetism at half-filling, when on-site repulsion 
terms of Hubbard-type are furthermore introduced. The 
model is expected to exhibit a quantum phase transition 
between the superconducting and the antiferromagnetic 
states, which is an essential feature of the cuprate super- 
conductors, when parameters are varied away from an 
exactly solvable point in a parameter space. 

Let us define the model. Let A be a rectangular lat- 
tice of the form A = [l,Li] x [1,^2] H 7? with periodic 
boundary conditions. It is assumed that Li is an odd 
positive integer and L2 = Li + 2. (We need these condi- 
tions to prove the uniqueness of the ground state as we 
will show later.) We denote by Cx,a (cj, the annihila- 
tion (creation) operator for an electron with spin a =|, I 
at site X. They satisfy the anticommutation relations 

{''2:,(TJ Cy^T-} = {Cx^ajCy^-r} = and {cj, Cy^T-} = 5x,yScr,T- 

We denote by $0 a state without electrons and by A^o the 
electron number. 

The hopping part of our Hamiltonian is given by 



Hi 



hop 



where tx,y = (1 



4A'^)i if a; = y, tx,y - -2Xt if \x - y\ = 1, tx,y = 2XH, 
if \x — y\ = V2, tx,y = if |a; — y| — 2, and zero 
otherwise 0- Here, it is assumed that t > and 
— l/4<A<l/4. In the wave space, it is represented as 
i?hop = Y.keKT.a=],i^i^)A,a^k,a whcrc e{k) = tg^{k) 
with g{k) = 1 — 2Acosfci — 2Acosfc2 for k — {ki,k2), 



(a) 





FIG. 1: (a) Lattice structure, (b) Dispersion relation. 
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The lattice structure and the single-electron dispersion 
relation are shown in Figs. H (a) and (b). 

Let us introduce new fermion operators corresponding 
to the single-electron states localized in the vicinity of 
site a; € A as follows: 



be,x,cr — 



y^A:\y — x\—l 



y6A;|i/-a;| = l 



e-ie-(y-)c. 



(2) 
(3) 



with 9 e {a,P,j} where a = (0,0), /3 = (0, tt) and 
7 ~ (7r/2,7r/2). The interaction discussed in this paper 
is attraction between electrons in these localized states. 
The interaction part of our Hamiltonian is given by 

Hint.,e = -WY, J2 ^Vx^-obe,x,-aai^„a^x,a (4) 

2;eA(T=T,i 

with Vl^ > 0. One easily finds that the summand in is 
bounded below by — 4(1+4A^)VF, which is attained by the 
states of the form aj. ^ ■ ■ ■ '^o- This indicates that 
Hmt,e describes attraction between two electrons with 
opposite spins. 

The whole Hamiltonian of our model is given by 



Hfi — H] 



hop 



Hint,e +ve 4i.aCo,cr (5) 



where t;^ = if ^ = a and vg > otherwise. The last 
term is added for a technical reason to show the unique- 
ness. 

To state our main result, we need to introduce fur- 
ther notation. Let G be the Gram matrix for the a- 
operator Q whose matrix elements are given by (G)x,y = 
{aj, g., ay^o-}. By a straightforward calculation, one finds 
that G is regular and that its inverse matrix is given 

by {G-^)x^y = l/|A|Efc6;cff"'(fc)e''-^"""^- Thus, it is 
possible to define dual operators of the a-operator as 
ax,a = Z]yeA(G'""^)a,Kaa,o-, which satisfy 

{'^•E.o-J '^J/,t} — {^'x,aT'^y,T\ ~ 6x,y^(j,T- (6) 

Since {aJ,^$o}KeA spans the single-electron Hilbert 
space, the feg-operators Q are expanded as 



&e,x,(T — / (U0)y^xO,y,o 



yGA 



(7) 



Here, the expansion coefficients (U0)y^x are given by 
{Ue)y,x = {al^„,be^x^c,}- One finds that {Ue)x..y = 
{Ue)y,x = {U'eTyx for 6 = a,l3 while {Ug)x.y = 
-iUg)y,x^mi^foTe^-f. 




(b) 



(c) 



FIG. 2: Wave-vector dependence of gg(~k) / g(k). (a) = a ■ 
(0,0). (b) e = /3 = (0,^). (c) e = 7 = (7r/2, V2). 



Using (Ue)x,y, let us define 



(8) 



which are the creation operators for electron-pairing 
states. The main result in this paper is as follows: 

Theorem. Suppose A ^ and consider Hg with W — 
t/A and fixed less than 2|A|. When is even, the 
ground state is unique, has zero energy, and is given by 



$0. 



(9) 



For odd Nc the ground state has positive energy. 

Before proceeding to the proof, we discuss the proper- 
ties of Hint,e and pairing states. 

By using the Fourier transforms of the c-operator, the 
fermion operators oj, g. , ^ ^ are expanded as 



ik-x 



e c 



fee/c 



I I keic 



(10) 
(11) 



where geik) — 2(cos(fci + 62) + cos(fc2 -t- 6*2)) for k = 
(fcijfe) and 6 — {61,92). One also finds from that 
4,, - l//REfce;c.rHfc)e"=-^4,.- Substitution of 
this and into = X^^eA '^I T^^ ^ i' ^bich follows 
from and ijHl, yields 



keic 



g{k) '^'^^T^-fea- 



(12) 



The precise expressions of ge{k) are given hy ga{k) = 
2 (cos fci -f cos k2), 3/3 (fc) = 2 (cos ki — cos ^2), and g^{k) = 
— 2(sinfci +sinfc2) (see Fig. These mean that and 
correspond to anisotropic spin-0 pairs while corre- 
sponds to an anisotropic spin-1 pair. 

We find from 1)10(1 and Hll|l that -ffint,e is expressed in 
the wave space as 

^mt,e = -ii ^fe,fe',<?4+g.T^fe'-9a^fe'a^fc,T (13) 

I I k,k',qelC 

Wlk',, = W(g{k + q)gg{k' ~ q)gg{k')g{k) 

+ge{k + q)g{k' ~q)g{k')gg{k)). (14) 
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One notices that our interaction Hamiltonian expressed 
as above contains scattering processes of electron pairs 
with non-zero total momentum. It should be also noted 
that for scattering processes with zero total momentum, 
only which are discussed in mean- field- type arguments, 
the amplitudes Wl _^ ,j = 2Wg[k + q)ge{k -f q)ge{k)g{k) 
become either positive or negative, depending on values 
of q and k. Nevertheless, the ground states are the su- 
perpositions of products of the electron pairs with zero 
total momentum. 

In the case oi 9 = 7, if we consider a Hamilto- 
nian obtained by replacing 7Jint,7 with H^^^ ^ — 

-WJ2xeAT.a=t,i^\,x,ab'f,x,aal^„aa:,a, which is inter- 
preted as attraction between electrons with the same 
spin, the following states become ground states for A 7^ 

and W = t/4: ^^m^,,,n^,^ = (C^^t)^ 

with even positive integers iVcj and N^^i less than 

|A|, where the pairing operators are given hy (^^ ^ — 

g. For iVeJ ^ 7Ve,i, *7.W.,T.JVea 

has a finite value of {Nc^i — Nc^i)/2, an eigenvalue of the 
third component of the total spin. This means that the 
coexistence of ferromagnetism and spin-1 pair condensa- 
tion is realized in the ground states of H!y. It is noted that 
the fully-polarized pairing states $^^Ar^^^o and $7,0,^04 
are stable for the on-site repulsion or the ferromagnetic 
interaction. These results may have some relevance to re- 
cently discovered materials exhibiting the superconduc- 
tivity as well as the ferromagnetism [3.ITol|. 

In the following, we shall prove the theorem for 6 — 
/3, 7. The case oi 6 = a can be proved in a similar but 
slightly simple way. 

Proof of Theorem for 9 — We first note that, 

by using the a-operator, TJhop is rewritten as i?hop = 
t^^^ij al „ax,a- Then, using this representation of 
i?hop as well as = t/A, we obtain 

Hg = wY^ al„he,x-ahl^^^_^ax,a 

a;GAcr=T,i. 

+ve ^ cJ,<tCo,<t- (15) 

Since all the operators in the right hand side are positive 
semidefinite, a state which is annihilated by these opera- 
tors is a ground state, having zero energy. We show that 
this is the case for '^e,N^ in 10. 

It follows from ®, Q and {hl^^J^ = that 

Ce^e,x,i«^.T- (16) 

Noting that (Ui3)x,y and {U^)x.y are symmetric and anti- 
symmetric, respectively, with respect to the exchange of 
X and y, we similarly obtain h\ ^ |flx,iCe = Ce^g x 
These relations imply that $e,JVc is a zero-energy state of 



the first term in the right hand side of H15|> . Furthermore, 
using co,<T = l/y|A|I]a,gA(l - 4A)"^a^,o- and 

E = (17) 

which follow from straightforward calculations, we find 
that co,(jCe ~ CgCo.o-- This together with the above result 
leads to Hg^g^N^ — 0. Therefore, ^e,N^ is a ground state 
of Hg. To see that ^e,N^ is actually a non-zero state, one 
rewrites (g as 

Cl= E (<T-4T)Ca (18) 

x£A\{0} 

by use of lfT7|l . Since each set of {(aj, ^ — al ^)<i>o}a:eA\{o} 

^'^d {&e,x,<T*o}a;GA\{o} is linearly independent [ij, ^e,N, 
is non-vanishing. 

The representation l|18() of motivates us to introduce 
the following lemma, from which the other statements in 
the theorem follow. 

Lemma. Suppose A ^ 0. Any zero-energy state of Hg 
with 9 = (3,^, W = t/A and less than 2|A| (where Ne 
is not fixed) is expanded as 

E '^^(n(«i.T-4T)) (nCi)*o (19) 

AcA\{0} \xeA ) \x^A I 

where the coefficients 4>a satisfy (j)A = 4>A' for any subsets 
A, A' such that \A\ = \A'\. 

This lemma implies that the ground state energy for 
odd Nc is positive. Suppose that there are two linearly 
independent zero-energy states for fixed even iVc. Since 
both of these states must satisfy the statement in the 
lemma, we find that the one is always represented by the 
other, which contradicts the assumption. Therefore, the 
ground state for fixed even N^. is unique. ■ 

Proof of Lemma. The parameter 9 is assumed to be P 
or 7 in this proof. Let us define oIq ^ = co,o- and ^ — 
ax,<T — ao,(T for x € A\{0} and also define b'g q ^ = co,o- 
and bg^^ — bg^xM for x £ A\{0}. These new operators 
satisfy the anticommutation relations 

{al,,ax,A^{b'g[o^.,b'g^,J = So,x (20) 

for X E A. Furthermore, each set of {a'^ ^^o}xeA 

and {bg ^ cr^o}xeA is linearly independent and spans 
the single-electron Hilbert space. Thus, the collection 

of states = (UxeAO-'lv) {Ux&Bb'e\x,-v) ^0 with 

subsets A and B such that \A\ + \B\ = forms a com- 
plete basis for the A^e-electron Hilbert space. Here, the 
spin index v is fixed to either | or |. 

Let $ be an arbitrary zero-energy state of Hg with 
W ^ t/A. We first expand $ in terms of the basis states 

*(A,s) as $ = Ea,bcA«^(^,-B)*(a,s) ^i*'^ coefficients 
4'{A,B)- To be a zero-energy state, $ must satisfy co,cr$ = 



4 



and bg ^ _^ax,cr^ = for cr =t , i and x G A. From the 
former condition and H20|l. we find that 4'{A,b) = if 
is contained in either A or B, or both. From the latter 
condition for xq G A\{0} with a =[ and {a^^^., aj,.o-} ~ 
Sx.y for x,y £ A\{0}, we obtain 

^ x[xo e A, xo ^ B]sgn[a;o; A, B] 

yl,BcA\{0} 

x0(A,s)*(A\{,„},Bu{.o}) =0, (21) 

where sgn[- • • ] is a sign factor coming from exchanges of 
the fermion operators, and x[ "event"] takes 1 if "event" 
is true and otherwise. Since all the terms in the left 
hand side are linearly independent, we find (j)[A,B) — 
if G A in addition to ^ B. This holds for any 
X e A\{0}, so that only the terms with A^B such 
that A <Z B (Z A\{0} can contribute to the expansion. 
Taking account of the above results, we rewrite $ as 

= T.a,B(iA;\a\>\b\^'(A,b)^\a,b) ^ith new coefficients 
'^'(A B)- Operating co,cr and b\ |«xo,T on $ in this form 
and repeating a similar argument to the above, we find 
that $ is expanded in terms of with A, B such 

that A = B c A\{0}. 

Any zero-energy state is thus written as 

EacA\{o} '^^^Ja.a) ^^lere (jjA = '/>(A,A)- Then, 
we again consider the zero-energy state condi- 
tion b\^^^axa,i'^ — ^ f^o^ ^0 G A\{0} and derive 
conditions on (j)A- Here, it is noted that bg^xo,a- 
is expanded as bg^xo^a = J2yeA\{o}i^s)y,xoa'y^^- 
From this and the anticommutation relation 

{b'lx,a^(^^o,<y} = {bl,x,a^^^oA = iU8)x,xo for X G A\{0}, 
we deduce 

Y X[y(^ A, y' £ A]sgn[y, y'; A] 

AcA\{0}y.y'eA\{0} 

xO<^^*(AuM,^\{.'})=0 (22) 

where F^"^, = (Ue)xo.y{Ue)y' ,xo, and sgn[- • • ] is a fermion 
sign factor. Let us choose a subset A which does 



not contain a nearest-neighbor site y of Xq but does 
contain next-nearest-neighbor site y' in the same di- 
rection. (The sites xq, y and y' are in the same 
axis.) For this set of sites, Fy°yt is non-zero. Then, 

by checking the coefficient of ^]^.,f i r we have 

J o {AU{y},A\{y'})^ 

{ssn[y,y';A](t)A+sgn[y',y;Ay,^y](f)A^,^J = where 
Ax^y is defined for x £ A and y ^ A hy Ax^y — 
iA\{x}) U {y}. Since sgn[y,y';A] = -sgn[y' ,y; Ay,^y], 
we obtain (j)A = 4>Ay,^y- 

Repeating the same argument for all x E A\{0}, we 
reach the conclusion that (pA = 4>A' whenever \A\ = \A'\, 
which completes the proof of the lemma. ■ 

Now let us consider the case of A = at half-filling 
with the inclusion of the on-site repulsion. Here, we 
furthermore assume that vg = Q for all 6 and that Li 
and L2 are even integers. In this case the Hamilto- 
nian becomes H^jj° = WY,^Y.a cl,a^e,x-abl,^_^Cx,a + 
^^x ^x \^x i'^^A^-^A with J7 > 0, which is still positive 
semidefinite. At half-filling, zero-energy states of on- 
site repulsion term are given by (IlajeA '^•r o-x) '^o with 
Ox =Tii- Then, considering the zero-energy condition 
for the first term in the Hamiltonian, we conclude that 
the ground states of -ffg are two-fold degenerate and 
given by the Neel states, exhibiting antiferromagnetism. 

One can readily find that the Hamiltonian Bg does not 
possess spin rotational symmetry. In the case of = a, 
however, we can construct an isotropic model with the 
ground state @ as follows. Let us define B[^^^ g = 

^J2xial^^bg^xA - al^bg^x,i)ibl,xA'^^A " ^I.xA^-^a) for 
6 — a, 0. A straightforward calculation yields that 

ibixA^-^A-bl^x.i^-^AKl = Clibl^xA'^^A - bl,x,iax,i) El: 
and, since H^^^^ g is positive semidefinite, © remains 
the ground state of Hg + -ff/^t g with W = i/4 for 
9 = a, /3. Furthermore, Hg + H[^^ g is isotropic since 

it commutes with S^^l = J^xi'^l^'^'-^'^ ~ and 
•^tot — l^x'^x I'^^d lis]- A construction of an isotropic 
model for the spin-1 pairing case and detailed investiga- 
tion of perturbed models of ours in both spin-0 and 1 
cases are left as an interesting future study. 
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